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ACTIGRAPHY  PROJECT 

WALTER  REED  ARMY  INSTITUTE  OF  RESEARCH 
Submitted  by 

Biomedical  Engineering  Center 
Purdue  School  of  Engineering  and  Technology,  IUPUI 


Introduction 

During  the  third  quarter  of  activity,  Purdue  School  of  Engineering  and  Technology  representatives 
continued  to  process  data  files  received  from  WEAIR,  participated  in  two  project  meetings,  one  at 
WRAIR  in  Washington,  DC,  and  one  in  Indianapolis,  and  worked  with  WRAIR  representatives  to 
refine  the  project  objective.  The  system  design  deliverable  was  completed  and  is  included  with  this 
report.  Some  support  activities  for  WRAIR’s  upcoming  “Rosetta  Stone”  study  were  also  carried  out. 
Following  are  summaries  of  each  of  these  activities.  Working  on  the  project  during  this  time  period 
were  Prof.  Russ  Eberhart,  Principal  Investigator,  and  Mr.  Xiaohui  Hu,  Research  Assistant. 

Project  Meetings 

Project  meetings  were  held  at  WRAIR  in  March  1999  and  in  Indianapolis  in  April  1999.  Methods  of 
data  analysis  were  discussed.  At  the  meeting  in  April,  it  was  decided  to  further  refine  the  work  done  in 
phase  I  to  include  a  software  tool  to  predict  PVT  performance  (reaction  times)  based  on 
polysomnography  data  (sleep  scoring).  The  system  design  information  presented  at  the  April  meeting 
was  deemed  sufficient  to  meet  the  deliverable  requirement,  for  which  an  extension  until  May  14, 1999 
has  been  requested. 

System  Design  Deliverable 

The  primary  goal  of  Phase  I  of  this  contract  is  to  develop  a  software  analysis  tool  that  classifies  a 
subject  as  asleep  or  awake  from  actigraph  data.  Neural  network  tools  that  achieve  a  testing  accuracy 
of  93-95  percent  have  been  developed  using  particle  swarm  optimization.  Included  with  this  quarterly 
report  is  a  diskette  with  the  trained  neural  network  weights  that  achieved  these  results  (files  r201.wts 
and  r202.wts)  and  the  particle  swarm  based  neural  network  executable  code  (file  psonn99.exe).  Also 
included  is  the  Visual  Basic  demonstration  code  for  particle  swarm  optimization  (file  swarm98a.zip). 
The  demonstration  code  can  be  used  to  optimize  several  benchmark  functions,  which  are  included. 

Attached  as  Appendix  A  is  a  printout  of  training  and  test  results  that  includes  results  for  the  two 
above-mentioned  weight  files.  It  can  be  seen  that  the  performance  on  test  data  (not  used  for  training) 
was  about  95.6%  correct  for  r201.wts  and  93.2%  for  r202.wts.  Attached  as  Appendix  B  is  a  listing  of 
the  Visual  Basic  source  code  for  the  particle  swarm  demonstration. 

Data  Preprocessing 

As  of  May  6,  1999,  we  have  received  65  actigraph  data  sets  on  CD.  A  listing  of  these  data  sets  is 
attached  as  Appendix  C.  Also  listed  in  Appendix  C  are  the  four  sleep  data  files  we  have  received. 


Each  sleep  data  file  starts  with  the  sleep  period  of  day  3  at  23:00.  Seven  days  of  sleep  of  a  specified 
length  (3,  5,  7,  or  9  hours)  follows,  followed  by  three  days  of  recovery  (9  hours  per  night). 

During  the  April  project  meeting,  the  preprocessing  of  the  actigraph  and  sleep  files  was  reviewed  in 
detail.  A  flowchart  describing  the  preprocessing  is  attached  as  Appendix  D.  During  the  discussion  of 
the  flowchart,  it  was  pointed  out  that  actigraph  data  values  are  for  the  previous  four  seconds,  where 
sleep  scores  are  for  the  next  30  seconds  (or  16  seconds  during  SLT).  It  was  decided  to  shift  all 
actigraph  data  4  seconds  “to  the  left”  (earlier)  so  that  it  is  consistent  with  the  sleep  scores. 

One  item  was  not  resolved.  When  we  process  the  data  into  60-second  blocks,  there  are  occurrences  of 
“1  0”  or  “0  1”  (sleep  wake  or  wake  sleep).  We  need  to  know  whether  to  score  these  one-minute 
periods  as  sleep  or  wake.  WRAIR  is  supposed  to  get  an  answer  to  us  very  soon. 

Data  Analysis 

During  the  April  project  meeting,  it  was  decided  that  we  would  focus  on  prediction  of  PVT  data  using 
the  sleep-wake  scoring  and  the  WRAIR  performance  model.  IV e  will  use  particle  swarm  optimization 
to  evolve  the  performance  model  parameters. 

WRAIR  will  provide  the  performance  model  as  soon  as  a  confidentiality  agreement  in  the  form  of  a 
CRADA  is  executed.  IUPUI  signed  the  CRADA  on  or  about  April  20, 1999,  and  forwarded  it  to 
WRAIR.  We  have  not  received  the  executed  CRADA  as  of  May  6, 1999.  Receipt  of  the  WRAIR 
performance  model  is  necessary  for  us  to  proceed. 

TeleActigraph-Related  Activities 

Activities  related  to  Precision  Control  Design’s  TeleActigraph  (TAG)  system  were  carried  out  during 
the  quarter  that  were  designed  to  support  WRAIR’ s  upcoming  “Rosetta  Stone”  study.  These  activities, 
which  were  coordinated  with  Precision  Control  Design  staff,  included. 

1)  Writing  and  compiling  a  program  to  pre-process  raw  data  as  downloaded  from  the  belt  unit  into  an 
ASCII  data  file  with  one  line  per  reading.  This  program  was  provided  to  WRAIR  in  February 

1999.  .  • 

2)  Writing  and  revising  a  Matlab  script  to  load  a  data  file  and  display  various  attributes  including 
power  spectral  density,  and  normalized  amplitude  versus  frequency.  This  script  was  provided  to 
WRAIR  in  February  1999. 

3)  Acquiring  resting  (awake)  data  from  five  additional  “normal”  subjects,  from  which  amplitude  and 
frequency  characteristics  of  typical  wrist  tremor  (with  unsupported  elbow)  can  be  derived.  We 
now  have  data  on  12  normal  subjects.  This  data  will  be  useful  for  the  Rosetta  Stone  study. 

4)  Becoming  more  familiar  with  TAG  System  operation  so  as  to  provide  observations  and  guidance 
to  WRAIR  on  system  operation  during  the  Rosetta  Stone  study.  Observations  and  guidance  (in 
addition  to  that  in  that  last  Quarterly  Report)  include: 

a)  Real-time  data  acquisition  (instead  of  off-line  downloading  of  data)  still  seems  to  be  a  very 
good  way  to  acquire  data  from  the  belt  unit,  at  least  during  sleep  periods.  We  have  ordered  a 
12-bit  National  digital-to-analog  PCMCIA  card  for  use  in  a  laptop  computer  that  is  running 
Lab  View  for  Windows.  A  LabView  script  has  been  written  that  permits  monitoring  of  data  as 
it  is  acquired  and  stored,  which  has  several  obvious  advantages. 

b)  We  have  installed  a  vibration  testing  system  featuring  a  VTS600  600-watt  amplifier  and  a 
Model  VG  100  vibrator  (Vibration  Test  Systems,  Aurora,  Ohio).  It  should  be  noted  that  this 
system  was  received  “on  loan”  from  Dr.  David  Burr  at  the  Indiana  University  Medical  School 
at  no  cost  to  this  project.  We  use  a  Techtronics  waveform  generator  or  a  pulse  generator  to 
drive  the  system,  and  monitor  results  with  a  digital  scope  and  a  PC  with  Labview  installed. 


TAG  units  have  been  preliminarily  characterized  and  calibrated,  and  benchmarked  against  the 
BMA-32  actigraph  units  that  will  be  used  concurrently  in  the  Rosetta  Stone  study.  Results  of 
this  work  were  provided  to  WRAIR  and  PCD  staff  members  at  the  April  project  meeting.  It 
appears  that  one  count  from  the  TAG  unit  in  Mode  9  corresponds  to  an  acceleration  of  about  2 
milli-gravities.  Further  work  will  be  documented  in  the  next  quarterly  report. 

Some  concern  exists  with  respect  to  the  fidelity  of  data  being  received  from  the  TAG  units. 
Driving  them  with  sine  waves  with  less  than  one  percent  harmonic  distortion  results  in  output 
data  that  is  not  as  nearly  a  sine  wave  as  might  be  desired.  Output  data  examples  obtained 
when  driving  the  TAG  with  a  6mm  peak-peak  sine  wave  at  1,  2,  and  3  Hz  are  attached  as 
Appendix  E.  This  matter  is  being  investigated  with  the  full  cooperation  of  PCD.  We  are  also 
examining  our  test  configuration.  Nothing  is  being  ruled  out  at  this  time. 
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Date:  1999-03-10  14:09:58 

Neural  Net  topology:  7  inputs,  14  Hiddens,  2  outputs 
Training  is  stoped:  14000  iterations  or  0.025000  errors 
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total 


0.890000  (  356/  400) 
0.902500  (  361/  400) 
0.876660  (  462/  527) 
0.899431  (  474/  527) 
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0.931689  (  491/  527) 


Date:  1999-03-10  14:09:55 

Neural  Net  topology:  7  inputs,  14  Hiddens,  2  outputs 
Training  is  stoped:  16000  iterations  or  0.025000  errors 


TrainSet:  157200.tr 
TestSet  :  157200. te 

classOO 

classOl 

Weights  :  r^laO.WTS 
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total 


0.917500  (  367/  400) 
0.937500  (  375/  400) 
0.920304  (  485/  527) 
0.944972  (  498/  527) 


0.920000  (  368/  400) 
0.912500  (  365/  400) 
0.907021  (  478/  527) 
0.933586  (  492/  527) 


frmSwarm  -  1 

i  'pa  y»4~ 4  ft *i  Qwflnn  VB5  software  by  Russ  Eberhart 

■  ££ S  the  open  statement  on  line  40  points  to  the  correct  drrectory 
»  May  18,  1998,  graph  iGbest,  iLbest 
»  May  29,  1998,  added  fO,  fl,  f2,  f3,  f5 
Option  Explicit 

Option  Base  1  .  .  . 

Dim  iPOPSIZE  As  Integer  '  Population  size  _ 

Dim  iDIMENSIONS  As  Integer  '  Number  of  dimensions 

Dim  iPopindex  As  Integer  '  Index  for  population 

Dim  iDimindex  As  Integer  '  Index  for  dimensions 

Dim  f INITWT  As  Single  '  Initial  inertia 

Dim  fMAXVEL  As  Single  '  Maximum  velocity  alio 

Dim  nMAXITER  As  Integer  '  Maximum  number  of  iterat 

Dim  niter  As  Integer  '  Number  of  iterations 

nim  fPr.<?n  As  Sinale  '  Position  for  each  particle 

SmSwPoM)  as  Single  '  Temporary  updated  position  lor  each  particle 

nim  fVein  As  Sinale  '  Velocity  for  each  particle 

Dim  f  DumVel  ( )  As  Single  '  Dummy  velocity  vector: ^Lrticle 

Dim  fBestPos ()  As  Single  ’  Best  previous  position  for  each  particle 

Dim  fMaxPos  As  Single  '  Dynamic  raa9® 

SrVaMfaf^oSle  ^  SionTrSr  value  calculated  by  evalO 

Dim  f  PbestVal  ( )  As  Double  '  Best  error  value  over  time  lor  each  particle 

Dim  fERRCUTOFF  As  Double  '  Error  value  at  which  system  stops 

Dim  iUSEBETTER  As  Integer  '  Usebetter  is  kind  of  a  momentum 

Dim  iBetterO  As  Integer  ’  This  gets  set  in  program 

nim  i rhpst  As  Integer  '  Index  for  global  best  particle 

Dim  iLbest  As  integer  •  Index  for  local  (neighborhood)  best  particle 

Dim  iLOCAL  As  Integer  1  Neighborhood  size  specified  m  run  file 

Dim  iHOODSIZE  As  Integer  1  Neighborhood  size  used  in  program  # 

Dim  iHOODINDEX  As  Integer  '  Neighborhood  index  (of f set  from  Partic  ® 

Dim  iNeighborO  As  Integer  1  Popindexes  of  neighbors  resulting  from  iLOCAL 
nim  ^Onifile  As  String  f  Filename  for  output 

Dim  iXcoord  As  Integer,  iYcoord  As  Integer  *  x  &  y  coords  for  plots 

Dim  iFNCTNO  As  Integer  '  Function  number  to  be  optimized  (1-fl,  etc.) 

Dim  fErrValDim  As  Double  ’  Error  value  for  a  single  dimension 

Const  conPI  =  3.14159265358979  ’Define  Pi 

Private  Sub  cmdSwarm_Click ( ) 

One>n  "c •\mvdocs\vb\psoinput.txt"  For  Input  As  #1 

Input  #1^  iPOPSIZE,  iFNCTNO,  iDIMENSIONS,  fERRCUTOFF,  fMAXVEL,  fMaxPos,  _  . 
nMAXITER,  iUSEBETTER,  f INITWT,  iLOCAL,  sOutfile 

'  Population  size,  function  no.  to  be  optimized,  no.  of  dimensions, 

'  error  cutoff  value,  maximum  velocity,  maximum  position  .  .  .  . 

t  (ranae),  maximum  no.  of  iterations,  use  momentum  (usebetter) ,  initial 
’  inertia  weight,  neighborhood  (0=global,  even  no.  for  local  neighborhood), 
1  filename  for  output 
Close  #1 


APPENDIX  B 


If  iFNCTNO  =  6  Or  iFNCTNO 


5  Then  iDIMENSIONS  =  2  'Set  #  of  dimensions  for  f5  &  f6  fnctns 


Re  Dim 
ReDim 
Re  Dim 
ReDim 
ReDim 
ReDim 
ReDim 
ReDim 
ReDim 


fPos ( iPOPSIZE,  iDIMENSIONS) 
fTempPos (iPOPSIZE,  iDIMENSIONS) 
fVel (iPOPSIZE,  iDIMENSIONS) 
fBestPos (iPOPSIZE,  iDIMENSIONS) 
f DumVel (iDIMENSIONS) 
fErrVal (iPOPSIZE) 
f PbestVal (iPOPSIZE) 
iBetter (iPOPSIZE) 

iNeighbor (-iPOPSIZE  To  iPOPSIZE) 


picSwarml . BackColor  =  vbBlack 

If  iLOCAL  >  0  Then  iHOODSIZE  =  Int (iLOCAL  /  2#)  *  2#  Else  iHOODSIZE 


=  iPOPSIZE 


r\d.IlU.UIt.L_L  LC  ,  ,  . _ 

«  Randomize  the  positions  and  velocities  for  entire  population 


frmS'warm  -  2 


For  iPopindex  =  1  To  iPOPSIZE 

For  iDimindex  =  1  To  iDIMENSIONS 

fPos (iPopindex,  iDimindex)  =  Rnd  *  fMaxPos  . 

fBestPos( iPopindex,  iDimindex)  =  fPos (iPopindex,  iDimindex) 
fVel( iPopindex,  iDimindex)  =  Rnd  *  fMAXVEL  . 

If  Rnd  >  0.5  Then  fPos (iPopindex,  iDimindex)  -  -fPos (iPopindex, 
If  Rnd  >  0.5  Then  fVel (iPopindex,  iDimindex)  =  -fVel (iPopindex, 

* f rmSwarm. Print  "Init.  Pos  &  Vel  ",  fPos (iPopindex,  iDimindex)  _ 
,  f™l  (iPopindex,  iDimindex)  1  for  diagnostic  purpose 

Next  iDimindex 
Next  iPopindex 


iDimindex) 

iDimindex) 


'  Main  swarm  loop  here 
For  niter  =  1  To  nMAXITER 


'  Update  inertia  weight;  linear  from  fINITWT  to 
flnerWt  =  ((fINITWT  -  0.4)  *  (nMAXITER  -  niter) 


0.4 

/  nMAXITER) 


+  0.4 


For  iPopindex  =  1  To  iPOPSIZE  'MAIN  main  loop  starts  here 

•Setup  dummy  velocity  vector  for  current  population  member 
For  iDimindex  =  1  To  iDIMENSIONS  ■  .  , 

fDumVel (iDimindex)  =  fVel (iPopindex,  iDimindex) 

Next  iDimindex 


iBetter (iPopindex)  =  0 


Set  to  0  unless  new  Pbest  achieved 


Select  Case  iFNCTNO 
Case  0 

Call  evalfO (iPopindex) 
Case  1 

Call  evalfl (iPopindex) 
Case  2 

Call  evalf2 (iPopindex) 
Case  3 

Call  evalf 3 (iPopindex) 
Case  5 

Call  evalf 5 (iPopindex) 
Case  6 

Call  evalf 6 (iPopindex) 
End  Select 


'  returns  fErrVal (iPopindex) 

'  evals  fO  spherical  function  error 

•  evals  fl  Rosenbrock  function  error 

•  evals  f2  Rastrigrin  function  error 

•  evals  f3  Griewank  function  error 
’  evals  f 5  He  function  error 

'  evals  f 6  Schaffer  function  error 


If  niter  =  1  Then  .  J  , 

fPbestVal (iPopindex)  =  fErrVal (iPopindex) 

iGbest  =  1 
End  If 

If  fErrVal (iPopindex)  <  fPbestVal (iPopindex)  Then  'If  new  Pbest 
fPbestVal (iPopindex)  =  fErrVal (iPopindex) 

For  iDimindex  =  1  To  iDIMENSIONS  '  Reset  Pbest  location  vector 
fBestPos (iPopindex,  iDimindex)  =  f Pos (iPopindex,  iDimindex) 
Next  iDimindex 


If  fPbestVal (iPopindex)  <  fPbestVal (iGbest)  Then 

•  color (iGbest)  =  DEFAULTCOLOR  need  to  dim  these  variables 
iGbest  =  iPopindex 
'  color (iGbest)  =  GBESTCOLOR 
End  If 


If  iUSEBETTER  =  1  Then  iBetter (iPopindex)  -  1 
Eikj  '  End  new  Pbest  loop 

Next  iPopindex  ’end  MAIN  main  loop  for  gold  gbest  only 


For  iPopindex  =  1  To  iPOPSIZE  'update  velocity,  position,  graph  position 

'  Does  neighborhood  calculation  of  iLbest 


If  iLOCAL  >  0  Then 

For  iHOODINDEX  =  0  To  iHOODSIZE 

iNeighbor (iHOODINDEX)  =  iPopindex  -  (iHOODSIZE  /  2#)  +  iHOODINDEX 
'for  iPopindex  =  l,goes  from  0  to  2  for  iHOODSIZE  of  2 
i  from  -1  to  3  for  iHOODSIZE  of  4 

'  Now  wrap  the  ends  of  the  array 

If  iNeighbor (iHOODINDEX)  <  1  Then  iNeighbor (iHOODINDEX)  =  iPOPSIZE  +  _ 
iNeighbor (iHOODINDEX) 

If  iNeighbor (iHOODINDEX)  >  iPOPSIZE  Then  iNeighbor (iHOODINDEX)  -  _ 
iNeighbor (iHOODINDEX)  -  iPOPSIZE 
'Start  with  iNeighbor (0)  as  iLbest  and  try  to  beat  it 
If  iHOODINDEX  =  0  Then  iLbest  =  iNeighbor (0) 

If  fPbestVal (iNeighbor (iHOODINDEX) )  <  fPbestVal (iLbest)  Then  _ 
iLbest  =  iNeighbor (iHOODINDEX) 

Next  iHOODINDEX 
End  If 


If  iLOCAL  =  0  Then  iLbest  =  iGbest 

'  Update  velocity  vector  for  one  particle  Russ  Reduced  version 

For  iDimindex  =  1  To  iDIMENSIONS  'flnerWt  below  . 

fVel (iPopindex,  iDimindex)  =  (0.5  +  (Rnd  /  2))  *  fVel ( iPopindex,  iDimindex) 
2#  *  Rnd  *  (fBestPos (iPopindex,  iDimindex)  -  fPos (iPopindex,  iDimindex) )  + 
2#  *  Rnd  *  (fBestPos (iLbest,  iDimindex)  -  fPos (iPopindex,  iDimindex)) 


If  fVel (iPopindex,  iDimindex)  >  fMAXVEL  Then 
fVel (iPopindex,  iDimindex)  =  fMAXVEL 
Elself  fVel (iPopindex,  iDimindex)  <  -fMAXVEL  Then 
fVel (iPopindex,  iDimindex)  =  -fMAXVEL 
End  If 

Next  iDimindex 

'If  it's  going  the  right  way,  keep  going 
If  iBetter (iPopindex)  =  1  Then 

For  iDimindex  =  1  To  iDIMENSIONS 

fVel (iPopindex,  iDimindex)  =  fDumVel (iDimindex) 

Next  iDimindex 
End  If 

Graphics  Loop:  Graphically  plot  updated  positions 

'Erase  old  position  for  two  dimensions  of  one  particle 
iXcoord  =  Int ( (picSwarml.ScaleWidth  /  2#)  *  (1  +  _ 

10  *  fPos (iPopindex,  1)  /  fMaxPos) ) 
iYcoord  =  Int  ( (picSwarml .  ScaleHeight  /  2#)  *  (!■*"_ 

10  *  fPos (iPopindex,  2)  /  fMaxPos)) 
picSwarml . PSet  (iXcoord,  iYcoord),  vbBlack 

For  iDimindex  =  1  To  iDIMENSIONS  'Define  new  positions  for  all  dimensions 
fPos (iPopindex,  iDimindex)  =  fPos (iPopindex,  iDimindex)  +  _ 
fVel (iPopindex,  iDimindex) 

Next  iDimindex 

'graph  new  position  for  two  dimensions  of  one  particle 
iXcoord  =  Int ( (picSwarml . ScaleWidth  /  2#)  *  (1  +  _ 

10  *  fPos (iPopindex,  1)  /  fMaxPos)) 

iYcoord  =  Int ( (picSwarml. ScaleHeight  /  2#)  *  (1  +  _ 

10  *  fPos (iPopindex,  2)  /  fMaxPos)) 

picSwarml . PSet  (iXcoord,  iYcoord),  vbRed 

'If  local  neighborhood,  graph  iLbest 
If  iLOCAL  >  0  Then 

iXcoord  =  Int ( (picSwarml . ScaleWidth  /  2#)  *  (1  +  _ 

10  *  fPos (iLbest,  1)  /  fMaxPos)) 
iYcoord  =  Int ( (picSwarml . ScaleHeight  /  2#)  *  (1  +  _ 

10  *  fPos (iLbest,  2)  /  fMaxPos)) 
picSwarml . PSet  (iXcoord,  iYcoord),  vbYellow 
End  If 


ErmSwarm  - '  4 


Next  iPopindex  'end  velocity,  position,  graph  loop 

'Graph  2  dimensions  of  iGbest 

iXcoord  =  Int(  (picSwarml . ScaleWidth  /  2#)  *  (1  +  _ 

10  *  fPos (iGbest,  1)  /  fMaxPos) ) 
iYcoord  =  Int ( (picSwarml . ScaleHeight  /  2#)  *  (1  +  _ 

10  *  fPos (iGbest,  2)  /  fMaxPos)) 
picSwarml. PSet  (iXcoord,  iYcoord),  vbWhite 


'Terminate  on  sufficiently  low  error 

If  (fPbestVal( iGbest)  <  fERRCUTOFF)  Or  niter  >=  nMAXITER  Then 

MsgBox  "Iter:  "  &  niter  &  "  Err:  "  &  fPbestVal (iGbest)  ',  vbYesNo 
Open  sOutfile  For  Append  As  #2  'psoutfile.txt 

Print  #2  "Iter-  "  &  niter,  "  Best  val:  "  &  fPbestVal (iGbest) ,  Now 
Print  #2,  "Neighd.  size  (0=global) :  "  &  iLOCAL,  "  Fnctn#  "  &  iFNCTt 
"  #Dims  "  &  iDIMENSIONS 

For  iDimindex  =  1  To  iDIMENSIONS  _  .... 


&  iFNCTNO, 


Print  #2,  "Dim. 
Next  iDimindex 
Close  #2 
End 

End  If 


&  iDimindex, 


'END  PROGRAM 


"  &  fPos (iGbest,  iDimindex) 


Next  niter 
End  Sub 


'  next  niter,  end  main  loop 


Sub  evalf 6 (iPopindex)  ’  evaluates  f6  function  error:  f Val (iPopindex) 

Dim  fNum  As  Double,  fDenom  As  Double,  fF6val  As  Double 


fNum  =  (Sin (Sqr ( (fPos {iPopindex,  1)  *  fPos  (iPopindex,  1)  +  _ 

fPos (iPopindex,  2)  *  fPos (iPopindex,  2)))))  A  2  -  0.5 
fDenom  =  (1#  +  0.001  *  (fPos (iPopindex,  1)  *  fPos (iPopindex,  U  +  _ 
fPos (iPopindex,  2)  *  fPos (iPopindex,  2)))  2 

f F6val  *  0.5  -  (fNum  /  fDenom) 
fErrVal (iPopindex)  =  1#  -  fF6val 
End  Sub 

Sub  evalf 0 (iPopindex)  'Evaluates  spherical  fO  function  error 

fErrVal (iPopindex)  =  0# 


For  iDimindex  =  1  To  iDIMENSIONS 

fErrValDim  —  ( fPos (iPopindex,  iDimindex))  2 
fErrVal (iPopindex)  =  fErrVal (iPopindex)  +  fErrValDim 
Next  iDimindex 


End  Sub 

Sub  evalf 1 (iPopindex)  '  Evaluates  Rosenbrock  fl  function  error 

fErrVal (iPopindex)  =  0# 

For  iDimindex  =  1  To  iDIMENSIONS  -1  ...  •  •  .j  \\  a 

fErrValDim  =  100#  *  (fPos (iPopindex,  iDimindex  +  1)  -  (fPos (iPopindex,  iDimindex)) 

2)  A  2  +  (fPos (iPopindex,  iDimindex)  -  1)  A  2 

fErrVal (iPopindex)  =  fErrVal (iPopindex)  +  fErrValDim 
Next  iDimindex 


End  Sub 

Sub  evalf 2  (iPopindex)  '  Evaluates  Rastngnn  f2  function  error 

fErrVal (iPopindex)  =  0# 


For  iDimindex  =  1  To  iDIMENSIONS 

fErrValDim  =  (fPos (iPopindex,  iDimindex))  A  2  +  10#  -  10#  *  Cos (2#  *  conPI 

opindex,  iDimindex))  „  , 

fErrVal (iPopindex)  =  fErrVal (iPopindex)  +  fErrValDim 

Next  iDimindex 
End  Sub 


fPos (iP 


frmSwarm  -  5 


Sub  evalf3 (iPopindex)  ’  Evaluates  Griewank  f3  function  error 

Dim  fFirstDim  As  Double,  fSecondDim  As  Double 
Dim  fFirstTot  As  Double,  fSecondTot  As  Double 
fErrVal (iPopindex)  =  0# 
fFirstTot  =  0# 
fSecondTot  =  1# 


For  iDimindex  =  1  To  iDIMENSIONS 

fFirstDim  =  (fPos (iPopindex,  iDimindex))  A  2 
fFirstTot  =  fFirstTot  +  fFirstDim 

fSecondDim  =  Cos (fPos( iPopindex,  iDimindex)  /  Sqr (iDimindex) ) 
fSecondTot  =  fSecondTot  *  fSecondDim 
Next  iDimindex 


fErrVal (iPopindex)  =  (0.00025  *  fFirstTot)  fSecondTot  + 
End  Sub 


1# 


Sub  evalf 5 (iPopindex) 
fErrVal (iPopindex) 
0.3  *  Cos (3#  *  conPI 

0.7 

End  Sub 


f5  function  error 

A  2  +  2#  *  (fPos (iPopindex,  2))  A  2  -  _ 

*  fPos (iPopindex,  1))  -  0.4  *  Cos(4#  *  conPI  *  fPos (iPopindex, 


'Evaluates  He 
=  (fPos (iPopindex,  1)) 


2)) 


Private  Sub  Form_Load() 
frmSwarm.  Show 
End  Sub 


APPENDIX  C 


File s  E : \Work\WRAIR\Data\actfile.txt  4/23/99,  12:34:28AM 


Actigraph  data  file 

308act.xls 
330act.xls 
333act.xls 
337act.xls 
351act.xls 
371act.xls 
3*.xls  17 

515act.xls 
518act.xls 
521act.xls 
543act.xls 
553act.xls 
556act.xls 
5* .xls  16 

711act.xls 
714act.xls 
725act.xls 
739act.xls 
757act*xls 
760act.xls 
7*. xls  16 

905act*xls 
928act.xls 
947act.xls 
962act*xls 
965act.xls 
968act.xls 
9*. xls  16 

total  65  files 


309act.xls 

331act.xls 

334act*xls 

349act.xls 

369act*xls 

372act.xls 


516act.xls 

519act.xls 

522act.xls 

544act.xls 

554act.xls 


712act *xls 
723act.xls 
726act.xls 
740act.xls 
758act.xls 


906act.xls 
929act.xls 
948act •xls 
963act.xls 
966act.xls 


310act*xls 

332act*xls 

335act.xls 

350act.xls 

370act.xls 


517act.xls 

520act.xls 

542act.xls 

545act.xls 

555act.xls 


713act.xls 

724act.xls 

738act.xls 

741act.xls 

759act.xls 


927act.xls 

946act.xls 

961act.xls 

964act.xls 

967act.xls 


PSG  data  file 

372sleep.xls 

•555sleep.xls 

711sleep.xls 

961sleep.xls 

total  4  files 


Page :  1 


Actigraph  output 


